Create a table summarizing data layers within each MTBS fire. Layers were assembled and reduced in Google Earth Engine, then the summary file exported to Google Drive. Downloaded the summary table locally from Google Drive, then ran this script to visualize and explore the data.

UPDATE: This script has been updated to filter for only ponderosa-dominated fire events (NOW COMMENTED OUT: with greater than 50% GEDI coverage). Additional metrics have been added, such as % USFS land, fire severity homogeneity, and burn severity amounts. Code hidden for ease of reading.

## [1] "C:/Users/tyler/Documents/Field Planning/CAREER"
## Reading layer `fires_stats_20210414' from data source 
##   `C:\Users\tyler\Documents\Field Planning\CAREER\data\fires_stats_20210414.geojson' 
##   using driver `GeoJSON'
## Simple feature collection with 37 features and 14 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -107.5838 ymin: 38.22499 xmax: -104.8615 ymax: 41.03432
## Geodetic CRS:  WGS 84
## # of fires in dataset:
## 37
Dominant forest type # of fire events
lodgepole 4
ponderosa 22
spruceFir 11

Note that a substantial part of the Hayden Pass fire is in a wilderness area, making it impossible for us to fly drones over that part of the fire. However, this part of the fire would likely be difficult to get equipment into, so we probably wouldn’t be flying there anyway due to the effort involved.

## {Generating outputs}
## PONDEROSA DOMINATED FIRES w/ GEDI >= 50% coverage
## Number of fires:
## 22
Fire Name Year Acres GEDI % Lodgepole % Ponderosa % Spruce Fir % USFS %
RED FEATHER 2017 1582 13.0 0.4 79.5 0.0 100.0
HAYDEN PASS 2016 17978 61.7 2.0 21.5 15.0 90.6
LOWER NORTH FORK FIRE 2012 3436 93.8 0.2 31.2 0.0 0.0
WALDO CANYON 2012 20112 58.6 0.3 47.2 0.2 78.1
HIGH PARK 2012 90769 46.8 13.9 49.6 1.1 51.3
SPRINGER 2012 1666 14.7 0.0 83.9 0.2 96.8
CRYSTALFIRE 2011 2855 95.3 2.7 54.2 0.0 41.5
DUCKETT 2011 4705 41.2 0.4 17.4 0.0 67.7
FOURMILE CANYON 2010 5865 91.1 1.2 67.0 0.0 4.5
PICNIC ROCK 2004 9014 88.2 0.0 15.9 0.0 34.9
OVERLAND 2003 3232 35.6 1.8 79.6 0.0 58.3
BIG ELK 2002 4347 98.5 19.3 47.9 0.0 92.3
SNAKING 2002 2080 91.4 0.2 65.4 2.7 69.3
HAYMAN 2002 129417 25.4 0.9 57.4 0.6 86.2
SCHOONOVER 2002 2821 0.0 0.0 73.2 0.0 87.7
BOBCAT 2000 9131 69.8 1.4 75.9 0.0 73.7
HIGH MEADOWS 2000 9607 0.0 0.0 75.0 0.0 55.4
BUFFALO CREEK 1996 9792 85.5 0.0 9.4 0.0 97.0
#6 1989 1547 93.2 0.0 6.7 0.0 53.1
SUNNYSIDE 1989 1416 47.8 0.0 27.5 0.0 54.6
UNNAMED 1989 1305 43.6 0.0 43.7 0.0 59.2
CANYON 1988 1092 47.5 0.0 40.8 0.0 62.4

Explore the fire data

Map

## Reading layer `fires_stats_20210414' from data source 
##   `C:\Users\tyler\Documents\Field Planning\CAREER\data\fires_stats_20210414.geojson' 
##   using driver `GeoJSON'
## Simple feature collection with 37 features and 14 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -107.5838 ymin: 38.22499 xmax: -104.8615 ymax: 41.03432
## Geodetic CRS:  WGS 84

Fire year histograms

Three histograms, one showing specific years of the fires, and two others showing decadal bins (either ’00-’10 etc or ’05-’15 etc)

Fire size histogram

Fire size as a function of year

With percent of fire inside USFS lands & GEDI coverage

Severity data

Severity proportions and homogeneity metrics for selected fires, based on MTBS-established severity classes

SUMMARY GRAPHS

ACTUAL HOMOGENEITY VALUES IF WANT TO SEE